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ABSTRACT 



o 
o 



We use population synthesis to construct distributions of gamma-ray bursts (GRBs) for differ- 
ent proposed progenitor models. We use a description of star formation that takes into account 
the evolution of metallicity with redshift and galaxy mass, the evolution of galaxy mass with 
redshift, and the star formation rate with galaxy mass and redshift. We compare predicted dis- 
tributions with redshift and metallicity to observations of GRB host galaxies and find that the 
the simple models cannot produce the observed distributions, but that current theoretical models 
can reproduce the observations within some constraints on the fraction of fallback black holes 
that produce GRBs. 
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Subject headings: Gamma Rays: Bursts, Supernovae: General 



Introduction 



With the advent of the SWIFT satellite the number of gamma-ray bursts (GRBs) with multiwavelength 
observations of afterglows and host galaxies is growing rapidly. As the data increases, it has the potential 
to allow us to probe deeper into the question of t he progenitors of t hese outbursts. On the surface, much 
of the data seems contradictory (for a review, see iFrver et al. I (|2007l) h recent data has confirmed that the 
explosions behind long-duratio n bursts produce supernovae, but we now know that some don't, or if they do, 
they produce dim supernovae ( Fvnbo et al. 20061 Berger et al.ll200(lk GRBs are associated wi th star forming 



regio ns, but they appear to occur in lower mass galaxies than their supernova counterparts (jFruchter et al 



2006); and although many bursts occu r at metallicities below l/10th solar, the average metallicity of bursts 
may well be closer to 1/3 to 1/2 solar (|Prochaskall2007l ) . 



All of this data may not be as contradictory as it seems. For instance, the idea that many bursts occur 
in low metallicity environments need not mean that bursts are exclusively produced by low- metallicity stars. 
If we consider that star formation is biased towards higher redshift and lower mass galaxies than a typical 
L* galaxy today, both of which imply lower metallicities, we might expect such an observational bias even if 
there is no metallicity dependence for GRB progenitors at all. As our sample of well-studied bursts grows, 
it will become easier to extract reality from observational bias. 

Theory has long made qualitative predictions about these bursts. Bursts are believed to come from black 
hole forming stars and the stars that form black holes tend to be the most massive stars. Thus, it is not 
surprising that GRBs and supernovae do not occur in exactly the same conditions. Since the strong stellar 
winds seen at high metallicites will prevent stars from collapsing directly to black holes, it has long been 
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suggested the rate of lo ng duration bursts will be higher at lower metallicities (jFrver. Wooslev fc Hartmann 
1999 ; Fryer et "al1l2007 ). Because the 56 Ni yield of a GRB explosion depends sensitively on whether the black 



hole is formed directly or by fallback in a weak initial explosion, theorists have also known that bursts would 
not necessarily alw ays be accompanied by bright supernova outbursts ( Fryer. Wooslev fc Hartmann 1999t 



Fryer et al.ll2006bl ). But to actually extract information about GRB progenitors with the current data, we 



must develop a more quantitative theory. 

In this paper we perform population synthesis calculations to determine if the common progenitor 
scenarios for long duration GRBs can be distinguished based upon the observational data. We combine 
the star formation rate as a function of redshift, evolution of cosmic metal density with redshift, the mass 
metallicity relation for galxies, and the specific star formation rate versus galaxy mass. Using a Monte 
Carlo technique, wc sample these distributions to produce the expected population of gamma-ray bursts and 
supernovae as a function of metallicity, redshift, and even host galaxy mass. Depending upon the progenitor, 
these distributions differ and we can use observations of bursts to constrain the progenitor model. 

SectionOdescribes the population synthesis Monte Carlo code used to create the simulated distributions, 
the descriptions of star formation and host galaxy properties used, and the GRB progenitor models. Section[3] 
compares our simulated distributions to observed GRB rates. Section 2] discusses the implications for the 
nature of GRB progenitors and our ability to constrain theoretical models. 



2. Calculations 

2.1. Initial Conditions 

We calculate the observed population of gamma-ray bursts by introducing a Monte-Carlo code that 
pulls from a sample of stars based on what we know of star formation. The fate of a massive star depends 
upon its metallicity, its initial mass, and its initial spin. To calculate the distribution of gamma-ray bursts, 
we first must assume a distribution for each of these quantities. 

Since the average metallicity of the universe decreases with increasing redshift, the metallicity of our 
entire distribution of GRBs clearly depends upon the star formation rate as a function of redshift. The mass- 
metallicity relation for galaxies also requires that we take into account the mass distribution of galaxies with 
redshift and the star formation rate as a function of galaxy mass. 

For our simple model, we use the functional form of the star formation rate per unit volume as a function 
of redshift (z) from Watanabe et al.(1999) 

J Alog(l + z) for z < z p 

togpiWi - | Mog ^ + Zp) _ B ^ z _ Zp) for z>Zp W 

A, B, z p are parameters that describe the star formation rate. The star formation rises from the initial burst 
of star formation and peaks at z p and then drops precipitously as we approach redshift of 0. Our standard 
model for the star formation history is based on a fit of this simple model to the current observations 
(including the metallicity evolution we use based on the work of Hopkins & Beacom 2006 - see below): 
(A, B,z p ) — (1.0,0.5,2.5). The star formation history remains uncertain, but its peak is now believed to 
lie somewhere between 2 and 3. We run 2 additional models probing the dependence of our results on this 
history: (A,B,z p ) = (1.25,0.5,2.0) and (1.0,0.1,2.5). To get the rate of star formation in a given redshift 
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bin, we must also correct for the volume filling factor: 

AV(z) = Az(1.0 - Vl-0 + z) 2 x (1.0 + z)- 3/2 . 



(2) 



For the mass-metallicity relation for galaxies we use the prescription derived bv lTremonti et al.l (|2004f ) 
from the Sloan Digital Sky Survey (SDSS) galaxy sample, which has the functional form 



12 + log(0/H) = -1.492 + 1.847 log M* - 0.08026(logM*) 2 



(3) 



Erb et al. ( 20061 ) show that the shape of the mass-metallicity relation remains essentially the same out 
to redshifts of z ~ 2, varying primarily in the normalization of the total metal content. We assume that the 
total metallicity scales along with the log(0/H) ratio. This is of course not strictly true, but it is the best 
observational approximation available. We assume that this functional form is valid to higher redshifts and 
shift the normalization of the distribution according to the t otal metal density of the un iverse as a function 
of redshift. For the total metal density we adopt the form in lHopkins fc Beacom (|2006h . 



p* = (a + bz)/[l + (z/c^hroMoyr^Mpc- 3 
p* = 63.7p z 

Pz(z) = I p z dz 



(4) 
(5) 

(6) 



where p* is the rate of change of the stellar mass density with redshift, p z is the r ate of change of th e 
metal density, and the coefficients a = 0.017, b = 0.13, c = 3.3 and d — 5.3 are defined in lCole et al. ( 2001 ). 



We use the galaxy mass function for z < 4 from the GOODS-MUSIC sample (jFontana et al 2006). This 



i s fit w ith a Shechter function with redshift dependent M*, and a* described as decribed in 
(120061) : 

/(M galaxy ) oc log[10 M - M *(*Y (z) e-™ M - M * (z) 

where M is the log of the galaxy mass, M*(z) — 11.16 + Q.YTz — 0.07z 2 , and a*(z) 
a redshift of 4, we set the mass function to the value at a redshift of 4. 



Fontana et al 



.18 



(7) 

0.082z. Above 



The most difficult part of the problem is a description for the star formation rate as a function of galaxy 
mass at different redshifts. The star formation in the local universe is dominated by dwarf galaxies, but by 
redshifts of 0.7 large spirals account for more than half. We produce a rough estimat e of th e specific star 
formation rate (M Q /M*) with respect to galaxy mass from Figure 13 of iPapovich et al.l (|2006l ). We fit a line 
to the points in the figures and take this as the specific star formation rate in three redshift bins centered 
on the redshifts of the three panels of their Figure 14: 



log(SFR/M gaUxy ) 



oc 
oc 
oc 



1.035ZogM galaxy ifz < 0.5 
0.8947ZogM galaxy if0.5 < z < 0.7 
0.44Zo#M ga i aX yifl.5 < z 



(8) 
(9) 
(10) 



where M ga i axy is the galaxy mass. These approximations are crude, but at least give us a zeroth order 
correction to the metallicity distribution of GRBs due to the mass-metallicity relation. 
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For our simple stellar mass distribution, we assume a single power law on the mass (M) of the star: 

f(M) = M- r dM, (11) 

where T is set to 2.7 for our baseline case. 

At this point, the uncertainties in our understanding of the evolution of angular momentum is too great 
for us to make any calculations of this quantity. For the time being, we assume that there is always sufficient 
angular momentum to make a GRB. This provides an upper limit on the GRB rate. 



2.2. Fates of Systems 

With the initial conditions outlined above, we produce a population of stars with given metallicities. 
To determine which is a GRB, we must then convolve this result with our knowledge of stellar evolution. 
Unfortunately, the final state of a star differs for different stellar evolution codes. The primary physical 
uncertainties in these codes include our understanding of convection in the stars (especially when it is 
coupled to nuclear burning), pulsational instabilities, and mass loss. Understanding rotation adds a new 
layer of uncertainties and we will not study this further in the paper. Binary stellar models are few and far 
between, generally focusing on specific objects. We will only make rough estimates of binary progenitors in 
this paper. 

As a test case, we first introduce a simple model that assumes all stars above a certain mass limit 
produce GRBs. This ignores the fact that mass loss is required to remove the hydrogen envelope, and too 
much mass loss prevents the formation of a black hole. We will use this simple model as a gauge to test 
our more realistic progenitor scenarios: basic single sta r prog enitors, binary star progenitors, and the more 



recent single-star models suggestied bv lYoon fc Langerl ([2005) 



For our basic single-star models, we take the Heger et al. (2003) paper using the strict definition of the 
plots where the metallicity is given as a log plot with "about solar" taken for exactly solar and the lower part 
of the plot taken for 10~ 6 solar. Using these definitions, we can predict the fates of a star as a function of 
metallicity and mass. Recall that our minimum requirements for a GRB are that the star must (i) collapse 
down to form a black hole and (ii) lose its hydrogen envelop^]. 

Our estimate of the effects of binaries is to assume that binary systems always remove the hydrogen 
envelope of the collapsing star. So our constraint for making GRBs is simply that the star must collapse to 
a black hole. Here we still use the Heger et al.(2003) distribution, but we can now extend the formation to 
include those stars that would normally still retain a hydrogen envelope. Note that specific binary scenarios 
may well produce distributions that are very different than the one used in our simple choice, but this 
provides us with a first guess of the distribution. 

We also include the more recent model suggested by Yoon & Langer (2005) that found that high rotation 
in stars could lead to extensive mixing that allows the near-co mplete burning of the hydrogen envelope. We 
base our models on the predictions in their most recent resultsl Yoon et al. 200tjh . 



Finally, remember that we can form black holes in 2 separate ways: a direct collapse to a black hole or 
the delayed collapse to a black hole via fallback after a weak supernova explosion. Some fraction of weak 



1 We assume that the star has sufficient angular momentum to form a GRB in all cases and hence, all our GRB rates are 
upper limits. We discuss the repurcussions of this assumption more fully in the conclusions. 
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supernova explosions will form very little 56 Ni (Fryer et al. 2006) and we will differentiate between these 
two black hole formation mechanisms. 



3. Results 

For most of our studies, we will use the classic single GRB rate as the basis for our intuition. Figure Q] 
shows the fraction of total outbursts per Az = 0.1 bin for 5 types of explosion: Type II supernovae, Weak 
Type II supernovae, Type Ib/c Supernovae, Weak Type Ib/c supernovae, and Direct-Collapse GRBs. The 
weak supernovae will produce black holes. For single stars, the weak Type Ib/c supernovae may also produce 
GRBs. The binary rate is very similar except that binaries can eject the hydrogen envelopes of the primary 
star, producing more Type Ib/c supernovae, and hence, more potential GRB systems through fallback. Note 
that the peak GRB rate from direct collapse happens at a higher redshift than core-collapse supernovae, but 
the weak supernovae are not so different than the redshift distribution of type II supernovae. Strong Type 
Ib/c supernovae have the lowest mean redshift. 

We expect GRBs to be made up of both direct-collapse and fallback black holes, but determining 
which type of black hole formation scenario dominates the distribution is difficult to determine. First and 
foremost, we do not know exactly which stars have weak supernova explosions versus the number that have 
no explosion ( Frverlll999 ) . Second, we know that fallback black holes will be less efficient than direct-collapse 



black holes at forming GRBs. This is because the accretion rate in a fallback black hole system is lower 
(possibly significantly lower) than a direct -collapse black hole sy stem. The lower accretion rate, the weaker 



the power available to drive an explosion. iPopham et al.l (| 1999T ) argued that the accretion rate would have 



to be above 0.01M Q y 1 to produce a GRB if it is powered by neutrino^. Recall also that the rates for 
GRBs here are upper limits (assuming all stars have enough angular momentum). 

Figure [2] shows the metallicity of GRBs and supernovae versus redshift for 6 separate simulations. 
Not e that without the fallback black holes, the metallicity distribution of single and binary stars under 



the Hcgcr et al. (2003) model predicts essentially all bur sts occur be l ow a metallicity of 0.1, very similar 



to the Yoon et all ( 2006 ) prediction. This is because the Heger et all |2003h models argue that above this 



metallicity winds are sufficiently strong that massive stars only form black holes through fallback. With 
fallback black holes, metallicities can be nearly as high as the supernova distribution (indeed, at high redshift, 
single stars predict metallicities higher than the supernova rate) . A simple explanation of any result that lies 
in between the pure direct-collapse and direct-collapse plus fallback systems is that the fraction of fallback 
black hole systems that form GRBs is smaller than those of direct-collapse black holes. Measuring this 
system could tell us about how fallback occurs in weak supernovae. 

Figure [3] shows the metallicity distribution for our binary models for two separate star-formation histo- 
ries. Although the flatter star formation history allows bursts and supernovae to occur at higher redshift, 
the other changes are fairly minimal. Only the highest redshift bursts will really help us constrain the star 
formation history. 

Figure H] shows the metallicity distribution of all explosions as a function of redshift. The total GRB 
distribution could be the summation of the direct-collapse and some fraction of weak Type Ib/c supernovae. 
Binary systems are similar, except that their distribution could include the direct collapse plus both types 



2 This constraint probably also exists for magnetically-driven outbursts at some level. 
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of weak supernovae. In any event, we expect the GRB population to have at least a slightly lower mean 
metallicity than supernovae. Included in Figure U is a plot assuming showing the combined distribution 
assuming 2% of all fallback black holes produce GRBs in the binary scenario. This produces a mean 
metallicity for GRBs just a little above l/10th solar. This mean metallicity could be significantly lower 
if direct-collapse GRBs dominate the GRB rate. It is no surprise that the mean mass of GRB-producing 
galaxies is also lower than normal supernovae (Fig. [5]). 

Figure [5] shows the number of binary GRBs including those from weak SN II and Ib/c(top), single 
star GRBs including weak SN Ib/c (middle), and all supernovae (bottom) per A = 0.1 bin in redshift and 
log metallicity (log^/z©)), normalized to the number in the peak bin. Contours are plotted for bins with 
20, 40, 60, and 80% of the peak number. On top of these contours, we show the current measu rements of 



meta llicity and redshift of GRBs associated damped Lya absorbers available in the literature (jProchaska 



2007). Triangles represent lower limits, circles metallicity estimates with error bars. The peaks of the 
predicted distributions are similar, but the GRBs have a much larger extent to low metallicities. This is 
greatest for the single star models. These models require the largest progenitor masses. As expected, the 
production of GRBs peaks between redshifts of 2 and 3, since that is where our star formation rate peaks. 

Figure [7] shows the number of binary GRBs including those from weak SN II and Ib/c(top), single star 
GRBs including weak SN Ib/c (middle), and all supernovae (bottom) per A = 0.1 bin in redshift and log of 
host galaxy stelar mass [log(M/M )], normalized to the number in the peak bin. Contours are plotted for 
bins with 20, 40, 60, and 80% of the peak number. While the populations have similar peaks, a larger fraction 
of GRBs occur in sub-L* galaxies. The higher specific star formation rate in low mass galaxies produces 
IMFs populated all the way up to the high mass stars that produce GRBs, while high star formation rate 
massive galaxies do not populate the extreme masses as well, producing a larger volume of relatively poorer 
stellar populations. This bias also explains the bias towards low metallicity in GRB hosts, as lower mass 
galaxies also have a lower metal content. 

Figure [8] shows contour plots similar to Figures [6] and [7] but for the "combined" scenario of Figured] By 
only including 2% of fallback black holes the peak of the distribution shifts to considerably lower metallicities 
and host galaxy masses. This distribution provides a reasonable approximation to the observed metallicity 
distribution. 



4. Conclusions 

In this paper, we present the theoretical distribution of GRBs as a function of redshift and of metallicity 
for the basic GRB models. We use current distributions of star formation, including the relative star 
formation rates of different-sized galaxies. In this manner, we can directly compare the predictions of the 
basic progenitor models to that of the growing set of observational data. 

From figure O which shows the same observational data overplotted on figure O we can immediately see 
that the small sample of GRBs is definitely at odds with the predictions of our simple models. For instance, 
the unmodified Yoon & Langer scenario cannot account for all GRBs. Their model is unable to produce 
any significant number of GRBs above metallicities of 0.1 solar. If the metallicity of GRBs without firm 
determinations lie even modestly above their lower limits, then at least half the sample has metallicities above 
0.1 solar. Similarly, direct black hole scenarios for both the binary and single star models are also biased 
towards low metallicity because mass loss at high metallicity is too extensive in the models. These models 
do have the advantage of producing a relatively flat metallicity ditribution for low to moderate redshifts, 
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and there is no obvious trend of metallicity with redshift in the data. 

Binary and single star models that include all fallback black holes produce too many GRBs at high 
metallicity by reducing the minimum mass for producing a GRB sufficiently that the GRB distribution to a 
large extent mimics the total supernova distribution. This also can be ruled out by the observations. 

We can make our models fit the data by assuming that the relative fraction of GRBs from fallback black 
holes is roughly equal to that of direct collapse. This can occur either by arguing that only a small fraction 
of fallback black holes have accretion rates that are high enough to produce GRBs or by moving down the 
critical mass that delineates fallback black hole formation from direct collapse black hole formation. Either 
of these possibilities could easily be true, and unfortunately, at this time, theory has not made any definitive 
predictions on cither of these. 

Quantitative estimates of host galaxy mass are not available for a significant number of GRBs. However, 
we can examine morphologies and see that most GRBs occur in small, actively star forming galaxies with 
peculiar optical morphologies. This is consistent with a bias towards low mass hosts. Again, the fallback 
scenarios produce too many GRBs at high host galaxy mass, but a reduced fraction of fallback black holes 
could easily be made to fit the data. 

There are several caveats to our current analysis. The first is that the current observed sample is still 
small and the observed distribution may change significantly when we have a larger statistical sample. Until 
the observational biases for the data we have can be understood and removed, it is dangerous to make any 
strong conclusions. In addition, we have made very simplistic models for the fates of massive stars. Binary 
progenitor models are much more complex than the simple models we have used, and until we have firm 
predictions from these models, it is difficult to make firm conclusions about these progenitors. Likewise, the 
most recent single-star models will evolve with time, and perhaps the inclusion of more complete physics will 
alter them to fit the data better. Fortunately, this framework can accomodate any progenitor prescription, 
so it can be used to test more sophisticated models as they become available. 
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Fig. 1. — Fraction of total outbursts per Az = 0.1 bin for 5 types of explosion: Type II supernova, Weak Type 
II supernova, Type Ib/c Supernova, Weak Type Ib/c supernova, Direct-Collapse GRB. The weak supernovae 
will produce black holes. Note that the GRB rate from direct collapse happens at a higher rcdshift than 
core-collapse supernovae, but the weak supernovae are not so different than the redshift distribution of type 
II supernovae. Strong Type Ib/c supernovae have the lowest mean redshift. 
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Fig. 2. — Metallicity of GRBs (solid circle and lines) and Supernovae (open circle s and dotted lines) versus 
redshift for 6 separate simulations: our simple model using just a mass limit, the lYoon et alj ( 2006T ) model 
prediction, Single and Binary models with and without fallback black holes included. The error bars cor- 
respond to the range enclosing 90% of all bursts. The solid line shows the metallicity of an L* galaxy as a 
reference point. Note that withou t the fallback black holes, the metallicity distribution of single and binary 
stars under the Heger et al. ( 2003T) model predicts essentially all bursts occur below a metallicity of 0.1, very 



similar to the lYoon et al.l (|2006f ) prediction. But with fallback black holes, metallicities can be nearly as 
high as the supernova distribution (indeed, at high redshift, single stars predict metallicities higher than the 
supernova rate) . 
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Fig. 3. — The metallicity distribution as a function of redshift for our direct-collapse binary GRBs and 
supernovae using two different star formation histories: A,B,z p from equation Q] are (1.25,2.0,0.5) and 
(1.0,2.5,0.1) for the decay and flat histories respectively. Within these narrow bounds of the star formation 
history, the changes in the history are modest. 
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Fig. 4. — Fraction of explosions per metallicity bin (0.1 in log space) as a function of metallicity for single 
stars. The total GRB distribution could be the summation of the direct-collapse and weak Type Ib/c 
supernovae. "Combined" is the sum of the direct collapse and 2% of the weak Type Ib/c rate. Binary 
systems are nearly identical, except that their distribution could include the direct collapse plus both types 
of weak supernovae. In any event, we expect the GRB population to have at least a slightly lower mean 
metallicity than supernovae. This mean metallicity could be significantly lower if direct-collapse GRBs 
dominate the GRB rate. 
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Fig. 5. — Fraction of explosions per galaxy mass bin (0.1 in log space) as a function of metallicity for 
single stars. The total GRB distribution could be the summation of the direct-collapse and weak Type 
Ib/c supernovae. Binary systems are nearly identical, except that their distribution could include the direct 
collapse plus both types of weak supernovae. In any event, it is clear that it could well be that GRBs are 
produced, preferentially, in smaller mass galaxies. 
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Fig. 6. — Number of binary GRBs including those from weak SN II and Ib/c (top), single star GRBs 
including weak SN Ib/c (middle), and all supernovae (bottom) per A = 0.1 bin in redshift and log metallicity 
(log(z/z0)), normalized to the number in the peak bin. Contours are plotted for bins with 20, 40, 60, and 
80% of the peak number. 
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Fig. 7. — Number of binary GRBs including those from weak SN II and Ib/c(top), single star GRBs including 
weak SN Ib/c (middle), and all supernovae (bottom) per A = 0.1 bin in redshift and log of host galaxy stelar 
mass (log(M/M©)), normalized to the number in the peak bin. Contours are plotted for bins with 20, 40, 
60, and 80% of the peak number. 
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Fig. 8. — Number of single GRBs including those from direct fallback and 2% of weak SN Ib/c per A = 0.1 
bin in redshift and log [M/H] (left) and log of host galaxy stelar mass (log(M/Mo)) (right), normalized to 
the number in the peak bin. Contours are plotted for bins with 20, 40, 60, and 80% of the peak number. 
The peaks are at much lower metallicity and host galaxy mass than for scenarios that include all fallback 
black holes. The observed metallicity distribution also matches better than for all fallback black holes or 
direct collapse or Yoon & Langer models alone. 
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Fig. 9. — Observ ed metallicity of GRB associated damped Lya absorbers versus redshift from the literature 
(|Prochaskall2007T ) ■ Triangles represent lower limits. 



